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Abstract 

We investigate the effect of a direct pomeron coupling to quarks on inclusive jet production in 
DIS and photoproduction. The direct pomeron coupling generates a point-like contribution to the 
diffractive part of the structure function F2, which is analysed on the basis of the latest HI and 
ZEUS data. Our model assumptions for the pomeron structure are consistent with the measured 
data. 



1. Introduction 

In diffractive production of hadronic final states in 
ep scattering, the proton stays intact or becomes a 
low mass state. Between the direction of the proton 
remnant, which goes down the beam pipe, and the 
produced hadronic system there is no colour flow, which 
allows of the possibility to observe large gaps in rapidity 
between these directions. The experiments HI |l| 
and ZEUS j^j at HERA have measured the portion of 
diffractive events to be « 10% of all events - not only 
in photoproduction (Q 2 < 0.01 GeV 2 ) but also in deep 
inelastic scattering (DIS) (Q 2 > 10 GeV 2 ). 

This paper is organized as follows. First, we describe 
our ansatz to analyse diffractive ep scattering. In section 
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3, we consider the pomeron structure function and the 
diffractive part of F2. We find consistency with the 
latest HERA data. Finally, we analyse jet cross sections 
in diffractive inclusive photoproduction and in DIS. 

2. Model for diffractive jet production 

There exist various phenomenological models || f|, 
||, H to describe the above mentioned diffractive 
nonperturbative QCD phenomena quantitatively. We 
follow a widely spread assumption, where the proton 
splits off a colourless object called pomeron (IP), 
which has the quantum numbers of the vacuum. 
Then, if factorization holds, the proton vertex can be 
parametrized by a IP— flux factor that depends 011 f = 
(p — p') 2 , the momentum transfer to the pomeron, and 
xp, the fraction of proton energy, that it carries away. 

In fact, this has been done in the past by various 
authors, who fixed their parameters with the help of 
pp scattering data. Inspired by Regge phenomenology, 
Berger et al. |7j found for the pomeron flux 
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with the Regge trajectory a = cxq + a't and the residue 
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function 0p v (t) = Pp p (0)e bot , where a = 1 + e, 



0.085, a' = 0.25 GeV" 



58.74 GeV 



and b = 4.0 GeV . 

This definition of the factorization in a pomeron flux 
factor differs by a factor of § from the definition of 
Donnachie and Landshoff ||] . In addition, they included 
the Dirac elastic form factor of the proton, which is given 

by 



For the hadron-like part of the unknown parton 
density functions of the pomeron, G^/p, we propose the 
ansatz 
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51 + r ' 

/3 G s/ ^ (/?) = /3 G^ /P 09) = |j9 G u/P 03) , 
0G g/P (f3) = 6/3(1- (3) j^- 

and vanishing charm contributions. Here, f3 = x/xp 
with Bjorkcn x. These functions obey the sum rule 
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with 5 = 3.24 GeV - . Equation || may be regarded as 
the most natural way to define a flux factor Q and will 
be used in the following. The same definition of the flux 
factor was used by Ingelman and Schlein Q . They made 
a different parametrization with two exponentials of the 
following form 

f P/p (t,xp) = i^Aexp^+Bexp' 3 '*) (4) 
2 KXp 

with the parameters k — 2.3 GeV 2 , A — 6.38, a = 
8 GeV" 2 , B = 0.424 and /3' = 3 GeV~ 2 . Here, the factor 
of 1/2 comes in because of normalization to one proton 
vertex. 

For our purposes, the momentum transfer t = (p 
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to the pomeron has to be integrated out, since the pomeron coupling to quarks g§ [ 



and are, in our case, defined for an input scale of 
Qq = 2.25 GeV 2 . Theoretical motivation on the basis 
of nonperturbative methods can be found in [fijj and || . 
The analysis of the authors suggests a behaviour of the 
diffractive parton distributions between (1 — /3)° and 
(1 - P) 2 for -» 1. 

The parameter r describes the unknown ratio of 
the gluon to quark content of the pomeron. To get 
a first insight into the structure of our model, we 
restrict the number of free parameters and choose a 
value r = 3, which represents simple gluon dominance in 
the pomeron. We carry out the usual DGLAP-evolution 
to get the right Q 2 dependence of these functions jl0| . 
Several groups considered the possibility of a direct 

A direct 



proton remnant is (still) not tagged. We get 
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with 



t 2 = -m 2 x 2 p/(l - xp) 



We emphasize that recently published HI data Q 
and ZEUS data || confirm this x~ n factorization. The 
actual values for the exponent n are 
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1.19 ±0.06 ±0.07 HI [1] 
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which are comparable to the exponent 2a(0) — 1 = 1.17 
from Regge analysis. However, since the pomeron might 
not be a real particle, there could be problems with the 
interpretation of (g) Q. 

Surely, the pomeron is in some sense only a generic 
object, that serves to parametrize a nonperturbative 
QCD effect. Although it is not considered as a physical 
particle that could be produced in the s— channel, we 
employ the concept of structure functions for it. 



pomeron coupling corresponds to a ^-function term in 
the pomeron structure function and produces a leading- 
twist behaviour in the pr spectrum. Our purpose is 
to find criteria that allow us to see a direct pomeron 
coupling in the data |l5| . 

As a consequence, similarly to 77 scattering, the 
jIP — > qq cross section also contributes to the pomeron 
structure function. Here, we assume a direct vector 
coupling of the pomeron to the quarks with coupling 
strength c. This is not really justified with respect to 
the C parity. However, the Q 2 dependence, that is 
~ logQ 2 at low x, is only weakly dependent on the spin 
structure. This can be seen, for instance, if one replaces 
the vector coupling by a scalar one. To remove the 
collincar singularity, we introduce the regulator quark 
masses m q and obtain for the point-like (pi) part 
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Figure 1. Generic diagram for the diffractive ep scattering 
process with a resolved photon 7, a resolved pomeron IP and the 
hard subprocess HS. 



3. The diffractive contribution to F2 

To leading order in a s , only the quark distributions of 
the pomeron enter into the deep-inelastic IP structure 
function Q 2 ), which is 



(11) 



= J> 2 /? [G ?/J >(/3,Q 2 ) + G ?/J ,(/3,Q 2 ) 

+2G p q l /p (p,Q 2 ) . 

The comparison with preliminary 1993 HI data [1] 
and ZEUS data [§ is shown in figure H and figure H. 
If factorization holds, which is favoured by the 
experiments for a large range of (3 and Q 2 values, 
the data points are proportional to Ff ((3,Q 2 ), i.e. 
F 2 (P,Q 2 ) = kFf(f3,Q 2 ) with a constant k that is 
determined by the _ZP-flux factor: 
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Figure 2. The pi distribution function for various Q 2 values 
and our choice of regulator quark masses. The right picture 
shows the e 2 weighted sum that enters into . In this plot x is 
the variable (3 in (lid). 
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The point-like contribution for the three quark masses 
m u — nid — 0.3 GeV, m s = 0.5 GeV and m c — 1.5 GeV 
is plotted in figure |[ The right picture shows the e 2 - 
weighted sum of the pi contributions from u, d, s and c 
quarks to F® . The dents in the curves are caused by the 
charm threshold. At fixed regulator mass, the maximum 
of the point-like contribution is shifted towards x = 1 
with increasing scale Q 2 . 

To get a constraint on the coupling c under the 
condition r — 3, we analyse the recently published 
data of F 2 D , the diffractive part of the proton structure 
function. 



dxp I dtfp/ p (x P ,t) 
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The absolute normalization due to k is very sensitive 
to the integration bounds, x T p n and Xp ax , which are 
taken from the respective experiments (see figure [| and 
figure ^). Unfortunately, no experimental errors on 
them have been published yet. A small reduction of the 
integration interval would improve our normalization to 
the data substantially. 

We concentrate therefore on the discussion of the 
shapes of the data in comparison to our model. We 
find that our model fits the shape of the data curves 
well. Especially, the Q 2 -evolution is fitted better for a 
combined ansatz, i.e., quarks in the pomeron and pi part 
(solid curves in figure ||), than for the DGLAP evolved 
quark distributions (short dashed line in figure |3|) or pi 
part alone. An alternative possibility to fit the data has 
been represented in ||. 

Finally, we fold the pomeron structure function ( |TT| ) 
with the pomeron flux factor in eq. (||) to get the 
diffractive part of the proton deep-inelastic structure 
function F®(x, Q 2 ). The relation is 



F 2 D (x,Q 2 ) 



dxp I dt fp/ p (xp,t)Ff (x/x P ,Q 2 ) 
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The upper bound xq, = 0.01 is an experimental choice. 
In contrast to the analysis of the pomeron structure 
function, the variable Bjorken x is now fixed (instead 
of (3). In figure |5| we compare with 1993 HI data (2(J. 
Again, we emphasize the consistency of the choice c = 1 
(r = 3) for the direct pomeron coupling in our model 
and the data. We use this value in the calculation of the 
diffractive jet cross sections. 
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0.0003 < x„ < 0.05 HI 



0.00063 < x„ < 0.01 ZEUS 




Figure 3. Comparison of our model predictions with 
preliminary 1993 HI data /J?j/. An important feature of the 
point-like contribution (long-dashed line) is the filling up of the 
quark distributions ( short-dashed line ) at higher (3 where they 
become less dominant. This results in flatter curves for the 
combined distribution (full line) and a better fit to the shape of 
the data. 



4. Jets in diffractive inclusive photoproduction 
and DIS 

The calculation of the differential jet cross section , , a .> 

of the process depicted in figure |l| is straightforward. 
Here, px and y are the transverse momentum and 
rapidity of one outgoing jet. In the case of 
photoproduction, we perform the evaluation in the 
ep laboratory system where the rapidity is positive 
for jets travelling in the proton direction. We have 
then the usual factorization of the photon flux factor 
at the electron vertex Jig ]. As is well known in 
photoproduction, the photon is resolved or couples 
directly to the final-state quarks. For the photon 
particle density functions, we take the parametrizations 
of GRV @. 

A more detailed discussion can be found in Jig ]. 
However, here we have included the improved quark 
distributions of the pomeron due to the performed Q 2 - 
evolution with Q = pt, the transverse momentum of 
the considered jet. Further, we use in our analysis more 
actual data of the run in 1993. 

The differential inclusive one-jet cross section is 



Figure 4. Comparison of our model predictions with 
preliminary ZEUS data ||/, analogous to figure^ 
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Figure 5. (x, Q 2 ) compared to preliminary 1993 HI 
fed/. The dashed lines represent only the contributions from 
quarks in the resolved pomeron, while for the solid curves the pi 
contribution with c = 1 is included. 



obtained by integrating out all kinematic variables over 
the allowed ranges without regard to the rapidity of 
the second jet, while for the two-jet cross section, we 
demand explicitly that the second jet does not enter the 
cone that is set up by the first jet around the outgoing 
proton direction. 

The results are shown in figure ^. Note the large 
rapidity gap in the forward direction between the 
diffractive and nondiffractive parts. In our analysis 
this is controlled by the x^ ax = 0.01 cut. This 
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Figure 6. The rapidity distribution of the a) one- and b) 
two-jet cross sections for fixed transverse momentum 
Pt = 5 Ge V in the ep laboratory system. Here, y is defined to be 
positive for jets travelling in the proton direction. For 
comparison, the nondiffractive cross section obtained with 
CTEQ parametrizations of the proton structure functions is also 
shown (solid line). Since the 1993 event rates (data points) of 
HI te/l/ and ZEUS fea/ are not normalized to the luminosity, we 
can compare only the shapes. 



Figure 7. The px spectra of the a) one-jet ancLb) two-jet cross 
sections for fixed rapidity y = 0. Like in figure la, the 
nondiffractive cross section obtained with CTEQ 
parametrizations of the proton structure functions is also shown 
(solid line). Since the 1993 event rates (data points) of HI Jjljj/ 
and ZEUS tell/ are not normalized to the luminosity, we can 
compare only the shapes. For higher p^ values the contribution 
from a direct pomeron quark coupling dominates the resolved 
part (c=0) and is favoured by the shape of the data. 



corresponds to a y max « 1.2 and can be compared 
to the experimental value ~ 1.5. The lower 

limit in the y-distribution is due to the cut on a; 7 that 
enters the improved Weizsacker Williams formula and is 
numerically determined by the experimental conditions. 

The need of a direct pomeron coupling becomes 
clear, if one inspects the slope of the pt spectra in 
figure ^. As expected, with the direct coupling, the pr 
spectrum does not fall off so strongly compared to the 
resolved pomeron contribution (c — 0). The point-like 
component of the structure function, however, does not 
play a significant role in this discussion, which can be 
understood by the gluon dominance in our model (r = 3 
in eqs. (|)). 

A reduction of r would increase the quark content 



in the pomeron due to the sum rule, eq. (|9|). But this 
would be accompanied by a reduction of the coupling c 
to satisfy the bounds coming from the analysis of the 
diffractive part of F2 in section 3. 

In DIS, the photon is always direct. The jet cross 
section has been calculated for our model in the j*p 
cms. Jets with positive rapidity are travelling now in 
the photon direction. With p, the four-momentum of 
the proton, k the four-momentum of the electron, and 
q, the momentum transfer to the photon, we define the 
usual kinematic variables s = (p + k) 2 , Q 2 = —q 2 7 
W 2 = (p + q) 2 , x = Q 2 /(2pq), y e = pq/pk, z = p T & jW . 
Finally, we set £ = x + p\,j{%)z(\ — z)s), which is 
the fraction of energy delivered from the proton to the 
subprocess (see figure 0). 
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Figure 8. T7ie ?/ distribution and pt spectrum of the one jet 
cross section in DIS, evaluated in the 7*p cms for fixed 
Pt = 2 GeV and y = 4, respectibly. We impose the experimental 
conditions Q 2 > 10 Ge V 2 and W 2 > 140 2 GeV 2 wAic/i 
correspond to the data points from ZEUS KiM. 



The inclusive one-jet cross section is then given by 
d 2 a 



dydp 7 



1 A-WVe fa 

dy e / ax- 

a+b Ja/y E 



dx dy e dy dpx 



,(14) 



where the kinematic bounds follow from the require- 
ments Q 2 > Q 2 min and W 2 > W^ in : a = Q 2 nm /s, 
b = max((2p T coshy) 2 , W^ in )/s. 

In the resolved pomeron case, we have 
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The contribution of the direct coupling is 
d 4 a 



dx dy e dy dp 2 , 
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The functions h u = Mh 
the Mandelstam variables s + i + 
subprocesses. For ~/q —> qg, we have 



hi) and hi depend on 
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Like in the photoproduction case, absolute experi- 
mental data for the rapidity distribution or pt spectrum 
are not yet available to us. In figure ||, we compare 
the shape of the pt spectrum with 1993 ZEUS data 
[ p4| . The experimental conditions were Q 2 > 10 GeV 2 , 
W > 140 GeV and 0.04 < y e < 0.95. No special case 
of our model is favoured, since the slopes of direct and 
resolved contribution are identical, but the shape can be 
approximately reproduced. 

In our analysis, we did not consider particle 
production or hadronization effects etc. The Monte- 
Carlo-programs POMPYT by Bruni and Ingelman ]r| 
or RAP GAP by Jung have been developed for a 
wide class of pomeron models and allow, together with 
other programs, the study of event characteristics. They 
are widely used by the HERA collaborations to interpret 
the large rapidity gap data. 

5. Conclusion 

In summary, we have studied the effect of a direct 
pomeron coupling on diffractive jet production at 
HERA. The concept of pomeron structure functions 
with DGLAP Q 2 -evolution has to be enlarged if the 
pomeron has a direct coupling to quarks. We have 
included the additional point-like part in the analysis 
of diffractive F2 data and find consistency for the 
assumption of a direct pomeron coupling to quarks. 
Some evidence for a direct coupling has been found 
in the pt spectrum of photoproduction. However, 
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our analysis is model dependent and second, except 
for the discussion of and , we have compared 
only the shapes and not the normalizations of the 
photoproduction and DIS cross sections. 
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